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Abstract 

In this paper, we give the parameter derivative representations in the form of 

„ , . n-l m k 


m=Q j=0 
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for some Koornwinder polynomials where A is a parameter and 0 < fc < n; n = 0,1, 2,... and present 
orthogonality properties of the parametric derivatives of these polynomials. 
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1 Introduction 

Recently, many authors [DElEllllEllIllilIHlIllIinilllllllllllllll] have studied the representations for the 
parameter derivatives of the classical orthogonal polynomials and various special functions which have many 
applications in applied mathematics, mathematical and theoretical physics and many branches of mathe¬ 
matics. In [aiiQ], the derivative of the Legendre function of the first kind, with respect to its degree v, 
[dP,^(z )G N), and its some representations have been examined by Szmytkowski, which are seen in 
some engineering and physical problems such as in the general theory of relativity and in solving some bound¬ 
ary value problems of potential theory, of electromagnetism and of heat conduction in solids. In explicit 
expressions of second-order derivative [d‘^Pv{z)/di3'^],j=o and of third-order derivative [d^Pu{z )/have 
been derived. In [miiiiii], the author has presented the derivatives of the associated Legendre function of 
the first kind with respect to its order and its degree and also a relationship between these derivatives. Such 
derivatives of the associated Legendre function are met in solutions of various problems of theoretical acous¬ 
tics, heat conduction and other branches of theoretical physics. In [allllElllllilIHlIIS], the representations 
of parametric derivatives in the form 

= ( 1 ) 

k=0 

for orthogonal polynomials in one variable, A being a parameter, have been studied. For instance, the 
representations of parametric derivatives have been obtained by Wulkow m for discrete Laguerre polyno¬ 
mials, by Froehlich [3] for Jacobi polynomials P^'^\x), by Koepf [4] for generalized Laguerre polynomials 
(x) and Gegenbauer polynomials (x) , by Koepf and Schmersau [S] for all the continuous and discrete 
classical orthogonal polynomials. In [5], Szmytkowski has derived again the expansions in the form of © 
for Jacobi polynomials, Gegenbauer polynomials and the generalized Laguerre polynomials by means of a 
method which is different from the methods given by Froehlich [5] and Koepf [1]. In [7], Ronveaux et al. 
have presented the recurrence relations for coefficients in the expansion 

=^^n,k{m,\)Pk{X',x) (m € N) 





which is more general than the expansion form of (HD- Moreover, Lewanowicz [6] has given a method to 
obtain iteratively explicit parameter derivative representations of order m = 1,2,... for almost all the clas¬ 
sical orthogonal polynomial families, i.e., continuous classical orthogonal polynomials, classical orthogonal 
polynomials of a discrete variable or q-classical orthogonal polynomials of the Hahn’s class. 

The classical Jacobi polynomials P^’^\x) are defined by the Rodrigues formula 




^(1 - x)-“(i +- xr+“(i+ 


and they satisfy the following orthogonality relation 

r pt^\x)pt^\x){i - xni+xfdx = 2 ^+11 

nlfa-I-p-I-2n-I-l)r(Q!-I-p-I-n-I-1) 


( 2 ) 


where 5n,m denotes Kronecker’s delta m- 

The generalized Laguerre polynomials defined by m 


hold 


= 


—rc n+Q; j 


dx^‘ 


[ L^^'>{x)L^“\x)e ^x°‘dx = + n + ^) ^ 

Jo n'. 

The representations of parametric derivatives obtained for the Jacobi polynomials (|3]) and generalized 

Laguerre polynomials (x) (0) are as follows 


dPlr’^\x) 
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n+k+a+B+l 




(/3 + l)n (2fc-|- Of -|-/3-|- 1 )( q ! -|- ,5 -|- l)fe p{a,p), X 

(a-I-/3-I-1)„ ^ (n - fc)(n-I-fc-I-a-I-/3-I-l)(/3-I-l)fc ^ 
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and 


dP^"'^\x) 

dp 


= E 
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n+k+a+P+1 




(g + l)n (—1)^^ (2fc + g +/3+ !)(« +/3 + l)fc p(Q,/3) / x 

(a-I-/3-I-1)„ ^ (n -/c)(n-I-fc-I-a-I-/3-I-l)(a-I-l)fc 


for a, /3 > —1 and 


dL!h\x) 1 r(a)/ X 


da 


k=0 


for a > — 1 where the Pochhammer symbol is defined by 


( 4 ) 


( 5 ) 


(a)o = 1, {a)k = a{a + 1)...(q; -I- fc - 1), fc = 1, 2,.... 

With motivation from the expansion (ED for orthogonal polynomials in one variable, we consider similar 
expansion in the form of 
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for orthogonal polynomials of variables x and y, with A being a parameter and 0 < k < n; n = 0,1, 2,.... In 
the recent papers [HE], parametric derivative representations in the form of ([6]) for Jacobi polynomials on the 
triangle and a family of orthogonal polynomials with two variables on the unit disc have been studied. The 
present paper is devoted to obtain parametric derivatives for the polynomials on the parabolic biangle, on 
the square and some new examples of Koornwinder polynomials introduced in m (see also [H]). Although 
the parameter derivatives of these polynomials with respect to their some parameters are in the form of (|6|). 
there exist some other parameters such that derivatives with respect to them are not in the form of (|6|). 
Because, some of the coefficients dn,j,m and depend on the variable x. 

The set up of this paper is summarized as follows. In section 2, we remind the method given by Koorn¬ 
winder m and some examples of this method. Section 3 contains parametric derivatives of Koornwinder 
polynomials on the parabolic biangle, orthogonal polynomials on the square, Laguerre- Jacobi Koornwinder 
polynomials and Laguerre-Laguerre Koornwinder polynomials. In section 4, some orthogonality relations for 
the derivatives of these polynomials are studied. 

2 Preliminaries 

First we recall some basic properties of orthogonal polynomials in two variables |20] . 

Let n be the set of all polynomials in two variables and let n„ denote the linear space of polynomials in two 
variables of total degree at most n. 

A polynomial p G n„ is called an orthogonal polynomial with respect to the weight function uj{x,y) if 

(p,?) := / P{x,y)q{x,y)u;{x,y)dxdy = 0 

J 

for all q G n„_i. Let Vn denote the space of orthogonal polynomials of degree n with respect to (,). 

In 1975, T. H. Koornwinder |19j constructed the following method to derive orthogonal polynomials in two 
variables from orthogonal polynomials in one variable. 

Let uJi(x) and u> 2 (y) be univariate weight functions defined on the intervals (a,b) and (c,d), respectively. 
Let p(x) be a positive function on (a, 6) which is either a polynomial of degree r, (r = 0,1,...) or the square 
root of a polynomial of degree 2r (r = i, 1, |,...). If p(x) is not a polynomial , c = —d < 0 and ^ 2 ( 1 /) is an 
even function on (—d, d). For fc > 0, let Pnix; k),n = 0,1, ... be orthogonal polynomial respect to the weight 
function and let qn{y),n > 0 be orthogonal polynomial with respect to the weight function 

012 ( 2 /) ■ Then, the family of polynomials 


Pn.k{x,y) = pn-k{x-,k)p’^{x)qk{—^), 0 <k<n 

p{x) 

are orthogonal with respect to the Koornwinder weight function (jj{x,y) = uji(x)uj 2 {-^^) over the domain 

fl = {{x,y) : a < X < b, cp{x) < y < dp{x)} 
with respect to the inner product 


/, 



fix, y)gix, y)oj{x, y)dxdy. 


Some examples of Koornwinder’s method are as follows: 

(i) Orthogonal polynomials on the parabolic biangle: 

For a, P > —1, Koornwinder polynomials on the parabolic biangle If = {(x,y) '■ y^ < x < 1} correspond 
with 


ujiix) = (1 — x)°‘x^, 0 < a; < 1, 

uj 2 {y) = {1 - y'^f, -l<y<l, 



p{x) = yj. 


These polynomials can be defined as 

y) = 0 < fc < n 

and they are orthogonal with respect to the weight function 

w(a;,2/) = (1 - xY{x - y^Y- 

In fact, 

{plYk\^^y)^Plnf\x,y)) ■■= J^Pn“f\x,y)P^f\x,y){l - xYix - y'^fdxdy 


= hYfSn^mdk,: 


where 




22/3+ir2(^ + /c + l)r(a + n-k + l)r(/3 + n + §) 


(n - k)\k\{2p + 2k + l)r(2^ + k + l)r(a + /3 + n + |)(a + /3 + 2n - A: + |) 


3\ ■ 


(ii) Orthogonal polynomials on the square: 

For a,l3,^,6 > —1, the polynomials defined by 


PlYk''^'^\^^y) = PlY-k\P)Pk'^\y), 0 < fc < n 


( 7 ) 


( 8 ) 


( 9 ) 


( 10 ) 


are orthogonal with respect to the weight function w(a;, y) = {1 — a;)“(l + xY{l — yY(\- + yY on the square 
O = {(x, y) : —1 < a; < 1, —1 < y < 1}. In fact, 


P 


(a,/3,7,5) / 


ix,y),P^f’'^'^\x,y)'^ := J^PY^f’'^'^\x,y)PYf’'^'^\x,y)il-xY{l + xf{l-yY{l + yYdxdy 

where dY’^'^ is given by ©• 


— Sl. ■ 

— ^n—k '^k ^n,m^k,j 


( 11 ) 


Some new examples of Koornwinder polynomials were introduced in m by using Koornwinder construction. 
These cases are as follows: 

(iii) Laguerre-Jacobi Koornwinder polynomials: 

The case of 

uJi{x) = x°'e~^, 0 < X < oo, 

uj2{y) = {1 - yY, -1 < y < 1, 

p{x) = X 

leads to the polynomials 


PlYfix, y) = LYYY''^^\x)x^PY’^\y-), 0 < fc < n 

which are orthogonal with respect to the weight function uj{x,y) = x°'~^e~^{x - 
domain Q = {{x,y) : —x < y < x,x > 0}. The following relation holds 


{piY^\x,y),PYf\x,y)'j : = J^P‘Yf\x,y)PYf{x,yY 


'(a 


— 5k - 

~ ^n,k kIn.mUk.j 


( 12 ) 

yY, (a,/3 > —1) over the 
- yYdxdy 



where 


( 13 ) 


(a,/ 3 ) 2 '^~*"^r(G; + Ti + /c + 2 ) 

^ (n-fc)!(^ + 2 fc + l) ■ 

(iv) Laguerre-Laguerre Koornwinder polynomials: 

From the Koornwinder construction with 

a;i(x) = x°‘e~^, 0 < a; < oo, a > —1 

W 2 (r/) = 0<y <oo, ^ > -1 

p{x) = X, a — (3 > —1, 

the Laguerre-Laguerre Koornwinder polynomials defined by 

0<k<n 

are orthogonal with respect to the weight function u}{x,y) = x°‘~^y ^over the domain 

n = {{x,y) : 0 < X < oo,0 < y < oo}. 

It follows that 

y): p!n,fi^^ y)) ■= Pn°‘f\x^ y)P^f\^^ y)x‘^-^y^e-^^+y/^^dxdy = t^^f^Sn,mSk,j 

where 

Aol,/3) r (/3 k l)r(a; + n + /p + 2) 


(14) 




k\{n — k)\ 


(15) 


3 Parametric derivatives of some Koornwinder polynomials 

In [HU], parametric derivative representations in the form of (jS]) for Jacobi polynomials on the triangle and 
a family of orthogonal polynomials with two variables on the unit disc have been studied. In this section, we 
derive parameter derivatives of Koornwinder polynomials on the parabolic biangle, on the square and some 
new examples of Koornwinder polynomials introduced in m (see also [H]). Since variable x is included in 
some of the coefficients, there exist some parameter derivatives such that they are not in the form (jS]). Now, 
we consider such representations of parameter derivatives. 

Theorem 1. For the Koornwinder polynomials over the parabolic biangle i^^y) defined by the 

parameter derivative with respect to the parameter a is as follows 


Ap(“’y) 

da 


{x,y) 


n — k — 1 
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p 




-t a + l3 + n + s + ^ 


{x,y) 


s=0 
n—k—1 


{a + fi + 2n — k — 2s — ^)(/3 -I- n — s -I- ^)s-i-i n(«,/3) / n 

“ (s-I-1 )(q! +/3-I-2n — A: — s-I-i)(a-I-/3-b n — s-I- 5 ) 5+1 


for n> k + 1, k >0 and 2 /) = 0 for n = k >0. 


(16) 






Proof. If we differentiate the both side of ([71) with respect to the parameter a, we get 


By using d^j), it concludes that forn>fc+l,A:>0 

n — k — 1 


d ^(a,P), . 1 


da 




a + /3 + n + s+| 

” ^ ^ {a + (3 + 2n — k — 2s — ^){fi + n — s + ^)s+i 


E 

s=0 


(s + l)(a + /3 + 2n—k — s + ^){a + fJ + n — s + 4)s+i 


Pi-P-i,kix,y)- 


It is obvious from d?]) that for n = fc > 0, ■^Pn%^\x,y) = 0. 

Theorem 2. For a, (3,^,5 > —1, the polynomials on the square defined by H1U\} satisfy 
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dP!rf’'^’^\x,y) 


da 


n — k + s + a + 13 + 1 




= i: 

s^O 

(a +/3 + 2n — 2fc — 2s — I)(/3 + n — fc — s)s+i p{a,p,'i,S), 

^ (s + l)(a + /3 + 2n-2fc-s)(Q; + /3 + n- fc - s)s+i ^ ’ 
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(17) 


dP. 


(a,/3,7.i5) 


n.k 


dp 


ix,y) ^ 

s=0 
n — k — 1 

+ E 

s=0 


for n > k + l,k > 0 and 


for n = k > 0. Also, 


n-fc + s + a + ^ + 1 ^ 

^ 2n - 2fc - 2s - 1 )(q + n - fc - s)g+i {a,p,^,S), , 

(s + 1 )(q! + /3 + 2n — 2fc — s)(q! + /3 + n — fc — s)s+i ri-s-i,k > V 


dP^°^n''^'^\x, y) dPk'^£''^’^\x, y) 


da 


dd 


= 0 




^7 


= E 

s=0 
k-1 


)(q;,/3,7,(5) 


7 + (5 + fc + s + l 


{x,y) 


(7 + 5 + 2fc — 2s — I)((5 + fc — s)s+l „(a,P,~i,S) , N 

+ ^ (. + l)h + < + 2t - .)h + i + t - .Wi 


s=0 


and 


9PnY’''’^\x,y) 


k-1 


dS 


= E 




7 + (5 + fc + s + I 


(a:,y) 


s=0 

fc -1 / ,Nfe_ 


, \ " (~1) "‘(7 + (5 + 2fc - 2s - 1)(7 + fc - s)s+l p(a,/ 3 , 7 ,^) 

^ (s + l )(7 + (5 + 2 fc-s )(7 + 5 + fc-s)«+i 


s=0 


for n > k > 1 and 


for n > 0 , fc = 0 . 


dP3:A''-‘\=:,v) dPl:fP-‘Hx,y) 


dj 


d6 


= 0 


(18) 


(19) 


( 20 ) 




















Proof. In view of equalities ([3]) and (H]), the proof is clear. □ 

Now, we can get similar results for Laguerre-Jacobi Koornwinder and Laguerre-Laguerre Koornwinder 
polynomials. 

Theorem 3. The representations of parameter derivatives with respect to the parameters a and fj for 
Laguerre-Jacobi Koornwinder polynomials defined by hl2\) are given by 


i: 


^—k—l 


P^°!’l^lix,y) 


g^-n,k ^ n-k-s 

s=0 

for n> k -\-l, k > 0 and -^Pn'^n^\x,y) = 0 for n = k > 0. Similarly, 

QP^n,k 

k-1 

+ E 


( 21 ) 




s=0 


(/3 + 2fc 2s 1)(^ 5 ) 8+1 s+l n(a+2s+2,/3) f ^ 


for n > k > 1 and 2/) = 0 n > 0,k = 0. It is seen that the parametric derivative with respect 

to the parameter /3 is not in the form of (0) since the coefficients include variable x. 

Theorem 4. For Laguerre-Laguerre Koornwinder polynomials defined by we have 

n—k—1 




s=0 


( 22 ) 


for n> k + 1, k >0 and ■^Pn'^n^\x, y) = 0 for n = k >0. Similarly, for n> k >1 


d 


fc-i 


1 




which is different from the form of (01 since the coefficients include variable x. Also, -§pP^Q^\x, y) = 0 for 
n > 0,k = 0. 


4 Orthogonality properties of parametric derivatives 


Now, we consider orthogonality properties for the parametric derivatives of the polynomials on the parabolic 
biangle, the polynomials on the square, Laguerre-Jacobi and Laguerre-Laguerre Koornwinder polynomials. 

Theorem 5. For Koornwinder polynomials on the parabolic biangle given by & and their derivative with 
respect to the parameter a, we have for n > k + l,k > 0, 0 < j < m; n, m G Nq 


p(a,P) d p(a,/ 3 )\ 

m,j ^ n,k / 


0 , 

0 , 


k = j,m > n 




n > m > k = j 


^n,k ’ 


k = j,n = m 


p(a,P) d p(a^0)\ 

m,3 ’ n,n / 


= 0 


and for n = k > 0 









where 


Aa,p) _ _ {a + j3 — k + 2m + |)(/3 + m+ |)n-m _t («,/?) 

n,k,m ^ ^ ^ ^ m,k 
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n,k 


n—k—1 
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1 

a + /3 + n + s + 


3 
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(«./3) 

n,k 


where is given by m- 

Proof. We will divide the proof into two cases. 

Case 1. We consider the case n = k > 0. It is seen that 


p(“./3) ^ p(a,0)\ 

m ,3 ^ n,n / 


= 0 


since ■§^Piifn‘\x, y) = 0 for n = fc > 0. 

Case 2. We assume that n > A: + 1, A: > 0. From (USD, we can write 


p(a,P) _^p(a,/3)\ ^ 


'i—k—1 


da 


Q. fj Tl S I 9 
s=0 ■* 


_ p(“./5) p(a./3) 
3 \'^m,j ’ ■^n,fc 


n—/c—1 

E 

s=0 


(g + /3 + 2n — fc — 2s — ^)(/3 + n — s + |)s+i /p(.a,p) p{a,i 3 ) 

(s + l)(a + ,5 + 2n - A: - s + 5 )(a +/3 + n - s + i)s+i \ 


(23) 


For this case, we consider three subcases. 

Case 2.1. Let consider the case k ^ j or k = j,m > n. It follows from (|SD that 

p(a,/3) p{a,/3)\ _ „ 

m,j ^ n,k / 


Case 2.2. Assume that k = j,n > m. Since the first inner product in the right-hand side of the equality 
(l23D from (jlD is zero, we get 

/p(o!,/3) d p(q,/ 3)\ _ ^ (g +/3 + 2n - A; - 2s - i)(/3 + n - s + i)s+i ,(a./3) c 

\^m,k ^ 0„^n,k /~ ifliO U 

\ uo. / (5 + 1)(q^ + p + ZTZ — k — s + 2 ) 1 ^ p n — s H- ^js+i 

which contains only one non-vanishing term with s = n — m — 1 for m > k. One may deduce that 
,(a,/3) __ {a + d — k + 2m -|- |)(/3 + m + ^)n-m _ 

n,k,m ^ ^ ^ ^ m,k 

where hpp is given by dSD. 

Case 2.3 Let k = j,m = n. Then, from the relation (jSD we have 

_p(a./3)\ _ 1 _/p(a./3) p{a,/3)\ _ _1_.(a,/3) 

\^n,k n,k / a + P + n + s + P ^ a + d + n + s+l, ’ 

which completes the proof. □ 

By using the parameter derivatives given in Theorem [5] and the relation dill) , the next theorem is readily 
verified. 










Theorem 6. For Koornwinder polynomials on the square defined by we have for 0 < j < m;n,m € Nq, 
n>fc + l,fc>0 


p(a,/3,7,(5) p(a,/3,7,(5) _ 

■mj ’ Qq, n,k ' 


0 , 


0 , 

k = j,m > n 

^(a,/3,7,5) 
^n,k,m ’ 

n > m > k = j 

^(a,/3,7,<5) 
-^n,k ’ 

k = j,n = m 


and for n = k > 0 


where 


p(a./3,7.<5) _^p(a,/3,7,(5) _ 

^m,j ’ ' 


0 , 

kf-j 

0 , 

k = j,m > n 

p(a,P,-y,S) 

^n,k,7n ’ 

n > m > k = j 

^(a,/3,7,5) 
^n,k ’ 

k = j,n = m 


/ p{a,P,l,S) p(a,/3,7,5)\ _ / p(a.^,7.<5) ^ p(a, 0 ,'f,S)\ _ q 

\ rn,j ) 5q; / \ / 


pi{(y.,0,S) _ “t“ 2??^ 2/c -t- l)(/3 k “t" ?TZ -j- j[a,l3) 7 ( 7 ,<5) 

(n - TO)(a +/3 + n + TO - 2fc + l)(a +/3 - fc + TO + l)„_m ^ 
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n —fe —1 

p,(a./3.7.5) ^ _ 

^ n-fc + s + a + /3 + l 

s=0 


j(«./3) j(7,'5) 

“n-fc “fe > 


p(a,/3,7,(5) _ 


(-l)”" + 5 + 2 to - 2 fc + l)(a - fc + TO + l)„_m Ja,p)j{j,s) 


(n — TO)(a + 5 + — 2fc + 1 )(q; + 5 — + l)ra-m ^ ^ 

where is defined as in m- 

Similarly, using the results in Theorem |3] and Theorem 01 one can easily obtain the next results. 

Theorem 7. For Laguerre-Jacobi Koornwinder polynomials P^^^\x,y) defined by \1‘A) . n,m £ No,0 < 
j < TO, we get for n>k + l,k>0 


p{a,P) d (a,p) 

rn,j 1 5q, n,fe 


G, 


0 , 

0 , 

n.k.m ’ 


k = j,m>n 
n > m > k = j 


and for n = k > 0 


where 


p{a,p) d p(a^p) \ _ „ 

rn,j ’ 5q, ' 


G 


»./3) 




n^k^m _ jy^^m,k 


where is given by il3\} . 










Theorem 8. For Laguerre-Laguerre Koornwinder polynomials n,m G No,0 < j < rn, the 

following results hold for n > k + l,k > 0 



'da 



0 , j 

0 , k = j,m > n 



and for n = k > 0 



where 


( 0 , 13 ) _ 



where t^jf'^ is defined by m- 


References 

[1] R. Aktas, A note on parameter derivatives of the Jacobi polynomials on the triangle, Applied Mathe¬ 
matics and Computation, 247 (2014), 368-372. 

[2] R. Aktas, On parameter derivatives of a family of polynomials in two variables. Applied Mathematics 
and Computation, 256 (2015), 769-777. 

[3] J. Froehlich, Parameter derivatives of the Jacoby polynomials and the gaussian hypergeometric function, 
Integral Transforms Spec. Funct., 2(4) (1994), 253-266. 

[4] W. Koepf, Identities for families of orthogonal polynomials and special functions. Integral Transforms 
Spec. Funct., 5(1-2) (1997), 69-102. 

[5] W. Koepf and D. Schmersau, Representations of orthogonal polynomials, J. Comput. Appl. Math., 90 
(1998), 57-94. 

[6] S. Lewanowicz, Representations for the parameter derivatives of the classical orthogonal polynomials , 
Rend. Circ. Mat. Palermo, Ser. II, Suppl. 68 (2002), 599-613. 

[7] A. Ronveaux, A. Zarzo, I. Area and E. Godoy, Classical orthogonal polynomials: dependence on pa¬ 
rameters,J.Comput. Appl. Math., 121 (2000), 95-112. 

[8] R. Szmytkowski, A note on parameter derivatives of classical orthogonal polynomials, arXiv: 
0901.2639v3. 

[9] R. Szmytkowski, On the derivative of the Legendre function of the first kind with respect to its degree, 
J. Phys. A, 39 (2006), 15147-15172 [corrigendum: J. Phys. A, 40 (2007), 7819-7820]. 

[10] R. Szmytkowski, Addendum to ‘On the derivative of the Legendre function of the first kind with respect 
to its degree’, J. Phys. A, 40 (2007), 14887-14891. 

[11] R. Szmytkowski, The parameter derivatives [d'^Pv{z )/and [d^Py{z)/dv^]y=o, where Pv{z) is the 
Legendre function of the first kind, preprint arXiv:1301.6586, 

[12] R. Szmytkowski, On the derivative of the associated Legendre function of the first kind of integer degree 
with respect to its order (with applications to the construction of the associated Legendre function of 
the second kind of integer degree and order), J. Math. Chem., 46 (2009), 231-260. 



[13] R. Szmytkowski, On the derivative of the associated Legendre function of the first kind of integer order 
with respect to its degree (with applications to the construction of the associated Legendre function of 
the second kind of integer degree and order), J. Math. Chem., 49 (2011), 1436-1477. 

[14] R. Szmytkowski, On parameter derivatives of the associated Legendre function of the first kind (with 
applications to the construction of the associated Legendre function of the second kind of integer degree 
and order), J. Math. Anal.AppL, 386 (2012), 332-342. 

[15] M. Wulkow, Numerical treatment of countable systems of ordinary differential equations, Konrad-Zuse- 
Zentrum Berlin, Techn. Rep. TR 90-8, 1990. 

[16] E.D. Rainville, Special Functions,The Macmillan Co., New York, 1960. 

[17] L. Fernandez, T. F. Perez and M. Pinar, On Koornwinder classical orthogonal polynomials in two 
variables, J. Comput. Appl. Math. 236 (2012), 3817-3826. 

[18] F. Marcellan, M. F. Marriaga, T. F. Perez and M. A. Pinar, Matrix Pearson equations satisfied by 
Koornwinder weights in two variables, arXiv:1411.2268 

[19] T. H. Koornwinder, Two-variable analogues of the classical orthogonal polynomials, in Theory and 
Application of Special Functions, R. Askey Editor, Academic Press, New York, 1975. 435-495. 

[20] C. F. Dunkl and Y. Xu, Orthogonal Polynomials of Several Variables. Encyclopedia of Mathematics 
and its Applications 81. Cambridge: Cambridge University Press, 2001. 


